Detection of the compressed primary stellar wind in 77 Carinae* 
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ABSTRACT 

A series of three HST/STIS spectroscopic mappings, spaced approximately 
one year apart, reveal three partial arcs in [Fe II] and [Ni II] emissions moving 
outward from 77 Carinae. We identify these arcs with the shell-like structures, seen 
in the 3D hydrodynamical simulations, formed by compression of the primary 
wind by the secondary wind during periastron passages. 

Subject headings: Stars: individual (t] Carinae) — stars: massive — binaries: general 
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1. Introduction 

rj Carinae (hereafter 77 Car) is now recognized to be a massive binary system in a highly 
eccentric orbit (e > 0.9). The primary star, 77 Car A, is in the luminous blue variable (LBV) 
phase, with a luminosity of 10 6 ' 7 l q and mass of about 120 M 0 (Davidson & Humphreys 
1997). It is surrounded by a massive stellar wind (M rs./ 10 3 M 0 yr 3 , Voo ~ 420 kms 1 ; 
Groh et al. 2012) that interacts with a companion, 77 Car B, that has not been detected 
directly yet. The physical parameters of 77 Car B are, thus, inferred from X-ray 
(M ~ 10 -5 M 0 yr -1 , ~ 3000 kms -1 ; Pittard & Corcoran 2002) and photoionization 

modeling (~ 10 5 L 0 , 34000 K< T eS < 3900 K; Verner et al. 2005; Teodoro et al. 2008; 
Mehner et al. 2010). 

The orbital orientation is such that the orbital plane has an inclination of i & 140°, 
with the semi-major axis pointing along position angle (P. A. 1 ) 314° and longitude of 
periastron oj ^ 240° (Gull et al. 2009, 2011; Madura et al. 2012). In this configuration, 

77 Car B is in front of 77 Car A at apastron, and the orbital angular momentum vector is 
parallel to the symmetry axis of the bipolar Homunculus nebula that surrounds the system. 

Projected on the sky, the secondary star directly illuminates and photo-ionizes the 
circumstcllar material within 1200 AU of the binary system primarily to the northwest, 
where the material is primarily approaching the observer. The companion star spends most 
of the orbit near apastron, so that the undisturbed wind from 77 Car A is located on the far 
side of the system, red-shifted as seen from Earth, and largely unaffected by the ionizing 
radiation from 77 Car B. 

During periastron passage, 77 Car B approaches within 1-2 AU of 77 Car A, leading to 
rapid changes in the shape of the wind-wind interaction surface and confinement of the 


X P. A. is defined as the angle, in the plane of the sky, eastward of North. 
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far-UV flux from 77 Car B. Hydrodynamical simulations (Madura et al. 2012, T. Madura et 
al. 2013, in preparation) show that, after each periastron passage, a highly-distorted volume 
of hot, low density secondary wind pushes outward into the slow, high-density primary 
wind, leading to the formation of a thin, high density wall surrounding the lower density, 
trapped wind of 77 Car B. This dense wall is accelerated to velocities somewhat higher than 
the terminal velocity of 77 Car A and expands both in the orbital plane and perpendicular 
to it, creating a thin, high-density sheet of trapped primary wind material. 

Using spectral imaging maps from the Space Telescope Imaging Spectrograph ( STIS ), 
on board the Hubble Space Telescope ( HST ), we identify structures which are consistent 
with these walls of compressed primary wind material. We identify three partial arcs formed 
by the close passage of 77 Car B around the primary star over the last 3 orbital cycles, and 
using the STIS maps, we determine their space velocity and age. 


2. Observations, data reduction, and analysis 

Observations using HST/STIS were accomplished through three guest observer 
programs, 12013 (2010 October 26), 12508 (2011 November 20), and 12750 (2012 October 
18), using the G430M and G750M gratings, centered at 4706A and 7283A, respectively, 
in combination with the 52" x 1" aperture. An extended region, covering 1" to 2" x 6.4", 
centered upon 77 Car, was mapped using 0.05" spacing offsets at P. A. constrained by the 
spacecraft. Exposures were approximately 30 seconds using CRSPLIT=2. An additional 7 
position sub-map, also centered on 77 Car, employed shorter exposures to allow for potential 
saturation of continuum on the central positions. 

Data were reduced with STIS GTO iDL-based 2 tools similar to the standard pipeline 


2 idl is a trademark of Exelis Visual Information Solutions, Inc. 
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reduction tools (Bostroem & Proffitt 2011), but with improved ability to spatially align 
the spectra to sub-pixel accuracy. For each emission line of interest, continuum levels 
were fitted for each spatial position along the aperture and a re-sampled data cube was 
generated, with the coordinates of right ascension and declination at 0.05" spacing, and 
Doppler velocity at 20 kms' 1 intervals. 

Since we are interested in structures formed in the primarywind around periastron, 
here we focus only on the following low-ionization lines: [Fe II] (Ao = 4729.39, 4776.05, 
4815.88, 7157. 13) 3 and [Ni II] (A 0 = 7413.65). Observations from program 12013 (2010 
October 26) did not include maps of [Fe II] A7157.13 and [Ni II] A7413.65. 

Visual examination of iso-velocity images in each data cube revealed multiple arcs in 
the low-ionization spectral lines used in this work. Sample images are presented in Figure 1 
for the [Fe II] A4815 transition. These arcs are visible between P. A. 90° and 200°. We 
identified at least three sets of partial arcs, labelled Al, A2, and A3 in Figure 2. The two 
innermost arcs, Al and A2, are conspicuously detected in all of the transitions from [Fe II] 
and [Ni II], The outermost arc, A3, is very faint. 

These arcs are not artifacts of the HST/STIS point spread function (PSF). The 
HST/STIS PSF diffractive rings vary slowly with wavelength, are symmetric about the 
central core and have amplitudes substantially weaker than the partial rings apparent in 
the images (Krist et al. 2011). As shown in the discussion below, these features correlate 
with multiple forbidden line structures seen at different wavelengths, expanding away from 
the central source. 

Comparison of the same iso-velocity image for different epochs (2010 October 26, 2011 
November 20, and 2012 October 18) revealed that Al, A2, and A3 move outward from the 


3 All rest wavelengths, Aq, are in A, in vacuum (Zethson et al. 2012). 
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central source (Figure 3). Because A3 is faint and hard to see in all the images, and Al is 
affected by uncertainties in the continuum subtraction method close to the central source, 
we concentrate our analysis on A2. We note, however, that all three arcs seem to expand 
at approximately the same rate, but some extended parts become very distorted or diluted, 
which makes the determination of accurate positions a difficult task. 

The position of the eastern portion of A2, relative to the central source, was measured 
using a radial flux profile cut along its direction (Figure 3). For each dataset (2010, 2011, 
and 2012), we extracted the flux profile between 0 " < r < 0.9" and 100° < P.A. < 105°, 
where the emission from A2 is strongest. 

We used three different methods for measuring the peak position in the radial flux 
profile curve. The first is a simple gaussian model, which is a good approximation when 
the flux profile of the shell is somewhat symmetric. The second method relies on the 
barycenter of the flux profile, and the third, developed by Blais & Rioux (1986), is based 
on the analysis of the derivative of the observed profile. The position of the peak is the 
median value of these methods. 

A comprehensive study comparing several methods to determine peak positions with 
sub-pixel accuracy, including the three methods used here, was done by Fisher & Naiclu 
(1996). The reader is referred to that paper for further details. Based upon the individual 
measurement errors, we estimate an accuracy of 0.25 pixel (0.013") for the position of 
the shell using the median value of the three methods mentioned before. The total error, 
however, may be larger than that because it is dominated by the dispersion in the measured 
position for various forbidden lines. 


3. Results 


We note that A2 is moving outwards with a proper motion /i = (0.05 ± 0.01)" yr -1 , on 
the plane of the sky (Figure 4), a value suspiciously close to the pixel size. However, many 
other regions across the held of view are fixed in position from image to image, giving us 
confidence that the images are all spatially registered correctly. Furthermore, the motion of 
portions of A2 depends on the local Doppler velocity (c/. Figure 4). Hence, the detected 
motion is not caused by a systematic misalignment of the images. 


The projected radius (r p ), at a given Doppler velocity (Vd), of a shell moving with 
space velocity V, during a time t, is given by 


r v {t) 


V t sin 


^arccos 



( 1 ) 


Taking the difference of the projected radius of A2, at two epochs separated by a time 
interval At = t 2 — H, leads to 


r p (t 2 ) - r-p(ti) = V At sin 



( 2 ) 


Setting the observations made in 2010 Oct. 26 as baseline, we have At = 1.07 and 
1.98 yr, for the two subsequent datasets. Thus, we used the measured projected position 
for A2 in the range —40 < Vp < +40 kms -1 (this is the small angle regime, where 
Vp/V -C 1) to estimate the space velocity, which resulted in V = (475 ± 29) kms -1 . The 
age of A2, t( A2), in 2010, is thus obtained from < r p (A2) > /V = (9.65 ± 0.55) yr, where 
< r p (A2) >= 0.41" is the average projected distance of A2, relative to the central source, 
at that epoch. 

Assuming that Al and A3 are moving at the same constant space velocity as A2, then, 
for a given epoch, we have 

r P (Al) = t(Al) 
r p (A2) “ £(A2) ’ 


( 3 ) 
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where r p (Al) is the average projected distance of Al from the central source, and i(Al) 
is its age. In 2010 and 2011, Al was too close to the central source. Thus, we used the 
observations of 2012 to measure the position and then estimate the age. In 2012, we 
have r p (Al) « 0.22", resulting in £(Al) = (5.98 ± 0.92) yr. The relation described in 
equation (3) also applies to A3, resulting in t( A3) = (17.13 ± 0.95) yr, since it was observed 
at r p (A3) ~ 0.6" in 2012. The rather large errors are dominated by uncertainties in the 
position due to the proximity to the central source for Al, and the weak emission of A3. 

Correcting for the time interval between the 2012 and 2010 observations, and comparing 
the ages in 2012, we finally have f(Al) = 5.98, £(A2) = 11.63, and £( A3) = 17.13 yr. This 
means that each shell is created approximately 5.6 yr after the earlier one, showing that 
they are directly tied to the orbital period of the binary system (P = 5.54 yr; Damineli et al. 
2008). 

To test for the consistency of our results across the range of Doppler velocities where 
A2 is observed, we adopted two simplified models in which we consider A2 to be a shell with 
either a spherical or an ellipsoidal geometry centered on r] Car, and allowed it to expand, at 
a constant space velocity of 475 kms -1 , during the corresponding age of A2 at each epoch 
of observations (£l exp , t2 exp , and f3 exp in Figure 4). 

At line of sight velocities near zero, both the spherical and the ellipsoidal models 
reproduce the observations very well, with standard deviations < 0.015". However, the 
spherical model does not match the observed position of A2 at red-shifted velocities 
(> 100 kms -1 ); A2 moves slower than predicted. We then considered the ellipsoidal model, 
to which we added some degree of asymmetry by changing the geometrical parameters (be., 
semi-major axis, eccentricity, orientation etc), improving the quality of the matching by a 
factor of 41% over the spherical model, for the full range of Doppler velocities. The final 
standard deviation for the ellipsoidal and spherical models is 0.034" and 0.048", respectively. 
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The ellipsoidal model better matches the observed positions because we can use two 
independent components for the space velocity: one along the ellipsoid’s major axis (fixed at 
483 km s” 1 ) and another along the minor one (300 kms” 1 ). It is important to remark that 
this is only a parametrization to reproduce the observed pattern of the Doppler velocity 
and is not a physical model of the geometry of A2. In our ellipsoidal model, when the angle 
of rotation ( 9 ) is zero 1 , the major axis is aligned in the plane of the sky; the minor is in the 
line of sight. We set the third axis to be the same size as the projection of the major axis 
onto the plane of the sky, the combination of which produces a circle projected on the sky. 
Therefore, our results suggest that A2 does not expand isotropically, i.e., with the same 
space velocity in all directions. 


4. Discussion and Conclusions 

We used the critical density for each forbidden line, calculated using the CHIANTI 
atomic database, to infer the typical local densities of the arcs. The outermost arc, A3, 
seen in the [Ni II] transitions but weakly in the [Fe II], must have a density in the range 
0.8 — 4 x 10 6 cm -3 . The strong [Fe II] emission from Al and A2, indicates that their 
densities must be > 4 x 10 6 cm -3 . 

Assuming that A2 is half a spherical shell (c/. Figure 5) with average radius of 0.4" 
(projected distance from the central source, at zero Doppler velocity, in 2010) and thickness 
of 0.05" (must be less than this since we cannot spatially resolve the shells), and using 
the critical density calculated before, we estimate a total mass of 2.6 x 10” 3 M 0 , for solar 
metallicity. 

Models for the primary star predict that its wind density should drop below 10 6 cm” 3 


4 The best overall match for the ellipsoidal model occurred for 6 = 8°. 
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for projected distances greater than 0.2" (Hillier et al. 2001; Groli et al. 2012). Since we 
detected arcs exceeding 10 6 cm -3 out to 0.6", this is naturally explained by the primary 
wind being compressed - a factor of 30 - by the secondary wind during periastron passage. 
This is also supported by the fact that the average age difference between the arcs is 
approximately 5.6 yr, nicely matching the binary system’s period. 

The primary mass-loss rate required to produce a shell with the same mass and 
geometry as A2, over a time interval equal to the binary’s period, is 9.3 x 10 -4 M 0 yr _1 , 
in agreement with the latest results from 2D radiative transfer spectral modeling 
(8.5 x 10” 4 M 0 yr _1 ; Groh et al. 2012). We note, however, that our result must not be 
taken as a strict lower or upper limit, since the actual geometry and thickness of A2 may 
be different from what we assumed. 

Simulations using 3D smoothed particle hydrodynamics (3D SPH) reveal compressed 
regions of the primary wind formed during the periastron passages and maintained 
for multiple cycles by the hot, low density, residual secondary wind (Figure 5). These 
simulations predict that these compressed regions accelerate from the primary wind terminal 
velocity (v^ = 420 kms -1 ; Groh et al. 2012) to about 460 — 470 kms -1 (T. Madura et al. 
2013, in prep.), which is consistent with our derived expansion velocity of 475 kms -1 for 
the shells. 

The departure from spherical symmetry for high red-shifted Doppler velocities suggests 
that the shells do not expand with the same velocity in all directions. This is explained 
by the fact that 3D SPH simulations show that material located in the orbital plane 
moves faster than material located off the plane. Thus, the observed Doppler velocity is 
a combination of velocities from different emitting regions along a specific viewing angle, 
which, in turn, depends upon the orientation of the shells on the sky, ultimately defined by 
the orbital orientation. 
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During periastron passages, SPH 3D simulations show that the secondary star creates 
a cavity deep within the primary wind, pushing a shell outwards. Since the low density 
secondary wind is moving at ~ 3000 km s' 1 , the interaction with the high density, relatively 
slow primary wind results in the formation of a high density, thin shell accelerated to 
velocities slightly higher than the terminal velocity of the primary wind. Hence, the 
geometry of the shell is strongly affected by the physical parameters of the binary system, 
and the projected appearance on the sky is defined by the inclination of the orbital plane. 

The latitudinal extent of the shells is controlled by the asymptotic opening angle, 9^, 
of the wind-wind collision region, given by (Usov 1992) 


and B refer to the primary and secondary star, respectively. 

Assuming fixed parameters for the secondary wind, a large primary mass loss rate 
decreases the bow shock opening angle 9^, so that the compressed primarys wind material 
is trapped near the orbital plane (Figure 5(c)); a low mass-loss rate increases 9^, and 
the compressed primary wind can extend well off the orbital plane (Figure 5(d)). This 
means that the geometry of the circumstellar arcs is an excellent diagnostic of the primary 
mass-loss rate: a lower mass loss rates produces a nearly continuous arc seen in projection 
around the star, while a higher mass loss rate produces a shorter, broken arc. 

The existence of a discontinuity in the arcs that we observe, more pronounced in Al as 
a lack of emission at P. A.« 135°, yields a lower limit to the mass-loss rate of the primary 
star. SPH simulations in 3D show that such broken rings are only apparent if the primary 
mass-loss rate is above 2.5 x 10 ~ 4 M 0 yr _1 , which suggests that, at least for the last event, 
the primary mass-loss rate must have exceeded this value. 



( 4 ) 


where 77 = Mg v^b/ Ma UooA is the ratio of the secondary to the primary wind momentum, 
M and Voo are the mass-loss rate and terminal velocity of the wind, and the subscript A 
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Our results also constraint the orbital orientation, with periastron occuring between 
200° < uj < 270°, i.e., during periastron, the primary star is between the observer and the 
secondary star. The opposite configuration, where the secondary, at periastron, is between 
the observer and the primary star, would produce blue-shifted shells seen to the northwest 
of the central source, rather than the red-shifted features to the southeast that we observe. 

Finally, we note that our results assume a constant shell space velocity. If these shells 
were accelerated at any time since their formation, their age, as derived here, would be 
overestimated, i.e., the shells would be younger than estimated, but the difference between 
the ages would remain the same 5.6 yr derived in this paper. 
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[Fe II] M815 (2012 Oct 18) 


-0.5 

(a) +160 km s 1 

(b) +200 km s' 1 

(c) +240 km s 1 

0.0 


*£■ 


0.5 



■ 

-0.5 

(d) +280 km s' 1 

(e) +320 km s -1 

(f) +360 km s _1 

o 

CD 




(arcs 

o 

o 

* 


ft 

to 

< 0.5 



■ 

-0.5 

(g) +400 km s -1 

(h) +440 km s' 1 

(i) +480 km s -1 

0.0 

ft 

n 

P, 

0.5 



1000 A.U. ■ 

1 — 1 


0.5 0.0 -0.5 0.5 0.0 -0.5 0.5 0.0 -0.5 

Aa (arcsec) 


~7 O ~7 a ~7 A ~7 O 


Fig. 1. — Iso-velocity images of [Fe II] A4815 for the observations of 2012 October 18. The 
velocities are indicated in each panel. The linear structures that appear in panels (g) and (h) 
are the result of a bad CCD column. The physical scale is based on an adopted distance of 
2350 pc (Smith 2006). We note that there are no extended structures at Doppler velocities 
higher than +440 kins -1 . 
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[Ni II] A74 13 (2011 Nov 20) 
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Fig. 2. — Iso-velocity image of [Ni II] A7413 at +200 kms -1 , obtained in 2011 November 
20. The arcs discussed in this work are indicated by the arrows and were labeled Al, A2, 
and A3. The dashed lines show the expected position of three spherical shells traveling at 
a constant speed of 475 kms -1 for 4.5, 10, and 15.5 years, for the inner, middle, and outer 
arcs, respectively. To determine the expansion rate, we restricted our measurements to the 
eastern component of A2. 



[Ni II] A.7413 (v=+200 km s' 1 ) 
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Fig. 3. — Flux profiles and iso-velocity images of two forbidden lines: (a) and (e) are for 
[Ni II] A7413 and [Fe II] A4815, respectively, along the position indicated by the arrow in 
panel (6), which has a length of approximately 0.9" and P.A. = 100°. The flux profiles 
were extracted from the iso-velocity images on the right. To illustrate the detected outward 
motion, in panels (c), (d), (/), (g), and (h) the curved solid line marks the position of A2 
in 2010 October 26. For better visualization, the flux profiles were shifted in the vertical 
direction. 
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Fig. 4. — Observed position of A2 as a function of Doppler velocity for the indicated 
epochs. Dashed lines show the projected position expected for a spherical shell expand- 
ing at 475 km s' 1 during £l exp , £2 exp , and t3 exp . The solid lines show the same, but for an 
ellipsoidal shell with space velocity of 483 kms -1 along the major axis, which is 8° off the 
plane of the sky, away from the observer. At zero Doppler velocity, both models have the 
same projected size on the sky. 
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Fig. 5. — 3D SPH simulations showing the formation and geometry of high density regions 
due to the compression of the primary wind by the secondary wind, in two different regimes 
for the primary mass-loss rate: 8.5 x 10” 4 M 0 yr _1 (left) and 2.5 x 10 -4 M 0 yr _1 (right). 
Upper panels: slices in the orbital plane (xy). Lower panels : slices perpendicular to the 
orbital plane (yz). In the high mass-loss rate regime, there is a gap in the polar region that 
is not observed in the lower mass-loss rate regime. The lower and upper limits of the density 
scale correspond to the critical density values for [Ni II] and [Fe II], respectively. 


